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Abstract. We analyse the quantum dynamics of a nanomechanical resonator coupled 
to a normal-state single-electron transistor (SET). Starting from a microscopic 
description of the system, we derive a master equation for the SET island charge 
and resonator which is valid in the limit of weak electro-mechanical coupling. Using 
this master equation we show that, apart from brief transients, the resonator always 
behaves like a damped harmonic oscillator with a shifted frequency and relaxes into 
a thermal-like steady state. Although the behaviour remains qualitatively the same, 
we find that the magnitude of the resonator damping rate and frequency shift depend 
very sensitively on the relative magnitudes of the resonator period and the electron 
tunnelling time. Maximum damping occurs when the electrical and mechanical time- 
scales are the same, but the frequency shift is greatest when the resonator moves much 
more slowly than the island charge. We then derive reduced master equations which 
describe just the resonator dynamics. By making slightly different approximations, we 
obtain two different reduced master equations for the resonator. Apart from minor 
differences, the two reduced master equations give rise to a consistent picture of the 
resonator dynamics which matches that obtained from the master equation including 
the SET island charge. 
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1. Introduction 

A very interesting class of nanoelectromechanical systems (NEMS) is composed of 
a mesoscopic conductor, such as a quantum dot, point contact, or single electron 
transistor, coupled electrostatically to a nanomechanical resonator [UEHHIIEIHIEIHIEI 
EJ EH HU E2J UH1 UH Um UIIl El- Such devices have potentially important applications 
as ultra-sensitive measuring devices PjJ HJ |H], as well as being interesting dynamical 
systems in their own right. A question of particular interest for these systems is under 
what circumstances the mechanical degrees of freedom require quantum mechanics for 
their proper description. 

When a resonator is coupled electrostatically to a mesoscopic conductor as a 
mechanically compliant voltage-gate, the position of the resonator affects the current 
flowing through the conductor and hence the latter can be used to monitor the motion 
of the resonator. However, the movement of electrons through the conductor necessarily 
acts back on the resonator affecting its dynamics in important ways. For weak electro- 
mechanical coupling, the back-action of the conductor on the resonator is typically 
analogous to an equilibrium thermal bath EH HH UE1 UE1 UH \H\ ■ The stochastic 
motion of electrons through the conductor gives rise to a force on the resonator which 
leads to fluctuations in its state that can be described by an effective temperature. 
Furthermore, the conductor can also damp the oscillations of the resonator, so that it 
reaches a thermal-like steady statej. 

The damping and heating effects of a conductor on a nearby nanomechanical 
resonator were first reported by Mozyrsky and Martin [3] who examined the quantum 
dynamics of a resonator coupled to a quantum point contact (QPC). In this device 
IT?)] the tunnelling matrix elements of a tunnel junction depend linearly on the resonator 
position. The back-action of the QPC on the resonator was obtained by tracing over all 
the electronic degrees of freedom to derive a reduced density matrix for the resonator. 
The resulting resonator master equation was found to be of the Caldeira-Leggett form 
and hence an effective temperature and damping constant arising from the back-action 
of the electrons could be ascribed to the resonator. 

The dynamics of a resonator coupled to a single electron transistor (SET) has also 
been studied theoretically [HI 13 EH E| • In the simplest case, electrons tunnel sequentially 
through a metallic island or quantum dot gated by a mechanical resonator [3 El- A 
classical master equation description of this system was proposed for this system in 
[7j and it was found that the resonator dynamics could be described by a Fokker- 
Planck equation ^Hj with the SET electrons again acting like a thermal bath. However, 
apart from weak electro-mechanical coupling, it was also assumed that the resonator 
moved very slowly on the time-scale of the electron tunnelling time and that the energy 
associated with the SET bias voltage was much greater than the resonator quanta. 
A subsequent quantum mechanical study of the closely related system of a resonator 

| Recently, it has been shown that more complex conductors like superconducting SETs can also give 
rise to negative damping of the resonator under certain circumstances |17| 
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coupled to a quantum dot also showed that the back-action of the electrons could 
give rise to resonator damping, even without assuming a slow resonator or very high 
bias voltage, but with a very different rate to that found in the classical study. 

In this paper, we begin from a microscopic description of the SET-resonator system 
and proceed to derive reduced master equations which describe the dynamics of the 
resonator alone. By first tracing over the microscopic electronic levels, we obtain a 
master equation which describes both the resonator and the excess charge on the SET 
island. We show that although this master equation is essentially equivalent to that 
proposed in the classical description jjj, it can nevertheless be used to investigate the 
interesting question of what happens when the resonator period matches the electronic 
tunnelling time. By solving equations of motion for the resonator moments, we find 
that the resonator is damped and reaches a thermal-like steady state even when the 
resonator period is of order or less than the electron tunnelling time, with maximum 
damping occuring when the electrical and mechanical time-scales are equal. 

Having obtained a master equation describing the SET island charge and resonator, 
we then go on to derive, via different approximations, two reduced master equations for 
the resonator alone. Although the two reduced master equations have differences in 
form, we find that (within their domains of validity) they both lead to a description of 
the resonator dynamics which closely matches that obtained from the master equation 
for the resonator and SET island charge. In this framework, the two different expressions 
for the resonator damping obtained before in quantum [Hj and classical treatments [7] 
emerge naturally as limiting values for the cases where the average tunnelling rate of 
the electrons is much less than the resonator frequency, and the opposite case where the 
resonator frequency is much smaller than the tunnel rate, respectively. 

The outline of this paper is as follows. In section |21 we outline the details of the 
nanomechanical-SET system we consider and derive a master equation which describes 
the dynamics of the resonator and the excess charge on the SET island for weak electro- 
mechanical coupling. In section |SJ we obtain equations of motion for the charge and 
resonator moments and hence explore how the resonator dynamics depends on important 
parameters such as the ratio of the relevant electrical and mechanical time-scales. Then, 
in section we trace over the SET island charge to derive a reduced master equation for 
the resonator alone, and compare the resulting resonator dynamics with that obtained 
in the previous section. Next, in section we use a rotating wave approximation to 
obtain an alternative reduced master equation for the resonator which is valid in the 
limit where the resonator frequency is much higher than the electron tunnelling rate. 
Again we compare the resonator dynamics predicted by this equation with the results of 
previous sections. Finally, in section H we present our conclusions and a brief discussion 
of our results. The appendices contain more details on certain aspects of the calculations. 
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Figure 1. Schematic circuit diagram for the SET resonator system. 



2. Master equation for resonator and SET island charge 



The nanomechanical resonator-SET system we consider is shown schematically in figure 
[D The SET consists of a metallic island, with total capacitance Cs, connected via tunnel 
junctions to two leads, with a bias voltage V applied across it. The nanomechanical 
resonator is adjacent to the SET island and is coated with a metal layer so that it 
forms a mechanically compliant voltage gate to which a voltage V g is applied. Motion of 
the resonator modulates the gate capacitance and hence the charging energy of the SET 
island, while changes in the charge on the SET island modulate the equilibrium position 
of the resonator. The resonator has a mass m and is treated as a single-mode harmonic 
oscillator with frequency u . Here we are primarily interested in the effect of the SET on 
the resonator and hence do not include any other influences on the resonator's dynamics 
in our description. 

Because the SET is metallic, the island as well as the source and drain leads contain 
many (microscopic) electron energy levels. The source (drain) lead has energy levels 
and a chemical potential fis(D) = —eV/2(+eV/2). The island has energy levels 6j and we 
set its chemical potential to zero for simplicity. The resistances of the tunnel junctions, 
Rs(D), are taken to be much larger than the quantum resistance, Rq = h/e 2 , so we can 
neglect higher order processes such as co-tunnelling [T%1 119j. 

The dynamics of the resonator is affected by the overall charge on the SET island, 
rather than the details of which of the microscopic electron levels are occupied, hence we 
introduce the macroscopic charge operator, N, for the total number of excess electrons 
on the SET island. Using this macroscopic charge variable, we can integrate over the 
microscopic electronic states, whilst still keeping track of the overall island charge |18j . 
The operator cf>, conjugate to N, can be used to form the operators e 1 ^ and e - ** which 
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increase and decrease the overall charge on the island respectively. 

The overall electrostatic energy of the SET island is determined not just by the 
island charge, but also by the gate voltage (which we write in dimensionless units as 
n g ) and the position of the resonator. We assume that the electro-mechanical coupling 
is weak, as has been the case in recent experiments jU IE] , and hence we consider only 
linear coupling between the SET island charge and the resonator position Bearing 
these details in mind, we write the Hamiltonian for this system as [T%1 ITTA !§]: 



where Eq = e 2 /2C*s is the charging energy of the SET island, a is the resonator 
annihilation operator, Ck, c q and q are the annihilation operators acting on the electron 
levels in the two leads and the island of the SET respectively, and Tu q )i are the tunnelling 
matrix elements between the microscopic states in the leads and the SET island. The 
coupling between the SET island charge and the resonator is given by x — rnu 2 xox q , 
where Xq is the shift in the equilibrium position of the resonator when an electron is 
added to the island (which in turn depends on details such as the gate capacitance 
and voltage as well as the resonator-island separation jjj), and x q = ^h/2muo is the 
zero-point position uncertainty of the resonator jTHj . 

Assuming that the charging energy of the SET island, Ec, is much greater than the 
thermal energy of the electrons, and that the drain-source voltage is not too large (i.e. 
eV/2 < Ec), then only the two consecutive charge states \Nq), \Ni) closest to n g will be 
involved in the dynamics of the system. In what follows we will restrict the charge state 
basis to just these two states, and introduce the convenient operators n = \Ni)(Ni\, 
a + = \Ni)(N \, and er~ = \N )(N 1 \. Choosing the origin of the resonator's position 
coordinate so that (a* + a) = corresponds to the equilibrium position when there are 
No electrons on the SET island, then the effective electrostatic energy of the SET can 
then be written (discarding constants) as Eq{1— 2n' g ), where n' g = n g —N [l — \ 2 hu /E c }- 
The parameter n' g ranges from to 1 and determines the difference in electrostatic energy 
of the two charge states. 

The form of the Hamiltonian (equation ^) is, apart from the leads, essentially that 
of the independent boson model 201 10| arid hence we make the canonical transformation 
usually applied to such systems to eliminate the term in the Hamiltonian describing the 
SET-resonator coupling. The operators in this canonically transformed picture are given 
by A = e - An ( at - a ) y 4 e An ( at - a ) ; where we have defined a dimensionless electromechanical 
coupling parameter, A = x/^o = xo/(2x q ). The canonically transformed Hamiltonian 
is then 




<-[ 



+ E C (N - n g ) 2 - xN{a ] + a) + hu a f a 
+ ^T fei (4qe-^ + c k c\e + ^) + ^T qi { 




(1) 



k q 



H = H s + H B 




k 
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cW"e A(at - a) + c q c\a + e- xiat - a) 



where the transformed system Hamiltonian is 
Hs = E c hn + hjjQO^a, 

with 



(2) 



(3) 



, ( 4 ) 

and the bath Hamiltonian describing the microscopic energy levels on the metallic leads 
and the island is given by 

(5) 



H 



B 



Converting to the interaction picture, Ar(i) = e^ Hs+HB ^ h Ae l ( H s+H B )t/h^ we obtain 

Hj = y^T u L f (t)Q(t)a-(t)e A ( at W- a W) +Q(t)4(t) ( x + (t)e- A(at( ' ) - aW) l ,(6) 



i 



where the sum over I runs over k and q. 

We wish to find an equation of motion for the density matrix of our system, p(t), 
which we define to be the resonator plus the macroscopic charge degree of freedom. For 
high resistance junctions, the tunnel couplings between the leads and the SET island 
(Tn) are weak, thus we assume that the microscopic energy levels in the leads and island 
remain in equilibrium throughout with occupancies set by the Fermi-Dirac distribution 
function (with the relevant chemical potential) so that = 2(0) for all time, t, 
where S(t) is the density matrix of the microscopic levels in the leads and the island 
(Born approximation [HIEJI22!)- We also make the conventional assumption that the 
full density matrix factorizes at t = 0, so that it takes the form p(0)S(0). Assuming 
that the electrons in the leads relax on a time-scale much faster than that over which 
pj evolves, we also make a Markov approximation [22] and hence arrive at a master 
equation of the form, 

t 

PM = - J? f dt ' TT [Hi(t'),Pim]} , (7) 



where the trace is over the microscopic electronic levels in the leads and island. 

Defining the bath operators B a = c»Cj and the system operator S 
can write the equation of motion for the system density matrix as: 

t 

1 



a - e A(at+a) we 



Plit) 



h 2 



dt>TrJ2\Tu\ 2 B u (t)S(t), Bl^^rpiit): 



1 

7^2 



dt'Ti \Tu\ 2 \Bl(t)S\t), [BafflSyifrit)! 



Li 



Notice that we need only sum over the microscopic levels once, as only matched pairs 
of creation and annihilation operators contribute when the electronic levels are in 
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equilibrium. The time dependence of the system operator S(t) = cr~(t)e A(at W~ a W) 
is somewhat complicated, and and so we write this as a sum of eigenoperators (defined 
by [Hs, S m ] = —Hu> m S m ), each with a time dependence of the form e _iWm< , where u m is 
an integer multiple of ujq. Anticipating the approximation we will make later that the 
electro- mechanical coupling is weak, we write out S(t) as a series in A: 

S(t)= ^a-(t)£ m e-^* 



a e 



1 + AaV" * 
A 2 



A 2 
2 



Xae- iuJot + — aVe 2 *"** 



A 2 A 2 A 2 

+— aae~ 2iuJot - —a)a - —aa j + 



(9) 

2 2 2 

The operators S m are products of different numbers of the operators a and a', with 
uj m determined by the number of each. Putting the time dependence explicitly into the 
equation of motion we obtain, 

t 



Pi®= -J-2 I dt'TrJ2\Tu\ 2 J2 e ~ 



Bug S v 



Bla+Slptt)! 



m,n 



s i(t— t')[ei-ei+E ch ]/H i(u m t— w n *') 



Bla+Sl, [Bua-S^p^ElllO) 



Li 



m,n 



We now assume that we can replace the integral over t' with an integral over r = (t — f) 
and that the resulting integrand is sufficiently peaked about r = that we can extend the 
upper limit to infinity (assumptions consistent with the Markov approximation |9*1 121j). 

oo 

The time integrals can be evaluated with the help of the expression: J dre lu)T = 

o 

7r5(u) + iPV (1/cjq). The principal value terms lead to coherent corrections to the 
evolution of the density matrix, but since we are working in the limit where the junction 
resistances are relatively high (Rs(d) Rq), these corrections are small and can be 
neglected CHI 123- 

Cyclicly permuting the bath operators under the trace, the master equation can be 
rewritten in the form: 



2tt 



■>ch 



l,i m,n 

+BuBkF [a + Sl(t),a-S m (t)]) pi(t)E 
where we have defined the super-operator, 



B\B a T [aS n (t),a+Sl(t)] 

(11) 



F\X,Y}p= l -{XpY + Y^pX^ 



- t {YXp + pX ] Y ] ). 



(12) 



We now perform the trace over the bath, using Tr[cjxfcS] = /(e^ — ps), etc., and 
convert the sums over the electron levels (i.e. over k, q, i) to integrals. After integrating 
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over 6^ and e q , the master equation becomes 

Pi{t) = 9d^2 / de if( e i)\ l - fi e i + E ch + fan - Pd)]F [ar~S n (t),ar + S} n (i)] pi{t) 

m,n 

+ 9dJ2 I de ^ 1 - f( e i)]f( e i + E ch + fan ~ Pd)T [v + Sl(t), a~S m (t)] p/(t) 

m,n 

+ 9s llJ de ^(^)[ 1 - f( £ i + E ch + ^ - Ps)\? [a-S n (t),a + Sl{t)] pj(t) 

m,n 

+ 9sJ2j deill-fie^fiei + E^ + ^-^la+Si^^-SUt^Piit), (13) 

m,n 

where the quantities gs{D) are defined by the relation 

9s{d) = p -J = -^iTfc^pDfe^A, (14) 

and it has been assumed that the tunnelling matrix elements Tu^ and the densities of 
states of the leads, D k ^, and island, Di are independent of energy over the relevant 
range of k(q),i. 

So long as the electron temperature in the leads is much less than the other 
relevant energy scales, the Fermi-functions can be replaced by simple step functions. 
Furthermore, for sufficiently large bias voltages (such that E D = eV/2 — E ch > tvui n and 
Eg = eV/2 + E c h > —fa n ) and sufficiently weak electro-mechanical coupling (we require 
A <C 1, so that the series can be truncated after a finite number of terms), tunnelling in 
the opposite direction to that set by the bias voltage cannot occur and hence the master 
equation can be written in the simplified form 

PM = 9s J2( E * + [°~ S n(t), ° + Sl{t)} p T (t) 



m,n 



g D Y,(E D - huj n )T [a + Sl(t), aS m (t)] p z (t). (15) 



111,71 



Converting the master equation out of the interaction and canonically transformed 
pictures, using a^ T = o~~ e A(a1 and oqt = a + Xn, and working to order A 2 , we 
find 

P(t) = - \ [Hg, p(t)] 

+ g s T [o-~(E s - \frw Q {a ] + a) + \ 2 huj ), a + ] p(t) 

-X 2 kj ),a-]p(t), (16) 

where 

H s = E c (l - 2n' g )n + hoj [a^a - An(a f + a)]. (17) 

The assumption that the electro-mechanical coupling is weak compared to the 
resonator quantum, A C 1 (i.e. x ^ fao) is the principal approximation made about 
the resonator in deriving this master equation. The derivation we have given here is 
convenient for later calculation, but it is also possible to derive equation (fTHjl under the 
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less restrictive weak-coupling approximation, x ^ e ^ \ as we discuss in Appendix A 



By casting the master equation in its Wigner function form (see Appendix A ), we see 
that it is essentially equivalent to the classical one proposed in [Zj, even for intermediate 
voltages [eV > fkj ) and fast oscillator motion (r <C 1). 



3. Resonator and charge dynamics 

The quantum master equation we have derived for the SET island charge and resonator, 
equation (jlfij) . is readily integrated numerically to give a complete description of the 
combined dynamics of the resonator and the macroscopic SET charge. However, 
for the normal state SET-resonator system we consider (unlike its superconducting 
counterpart [T7j) it is also possible to obtain a closed set of equations of motion for the 
moments of the electrical and mechanical degrees of freedom§. In this section we use 
a combination of analytical and numerical methods to investigate the dynamical and 
steady state properties of the resonator and SET moments. In later sections, we make 
further approximations in order to derive more compact analytical descriptions of the 
resonator dynamics and comparison with the results obtained in this section play an 
important role in judging their fidelity. 

We can calculate the equations of motion for the moments by multiplying the 
master equation by the appropriate operator and tracing over the resonator and charge 
states. Assuming for simplicity that gs = 9d = g, & closed set of equations for the first 
moments is obtained, 

j t (n) =g [eV/2 - E c (l - 2n' g ) + Xhu (a^ + a}} - geV(n) (18) 

—(a) + a) = iuj^io) — a) (19) 
at 

-l( t_ a \ =i Uo [(ot + a) -2A(n)l . (20) 
at 

These equations of motion are naturally equivalent to those for a classical system 
(ref. jZj). Here we investigate their behaviour over a wide range of parameter values 
and find in each case that the resonator motion is closely analogous to that of a damped 
harmonic oscillator with a shifted frequency. Of particular interest is the effect of varying 
the ratio of electrical and mechanical time-scales, V = geV/u>o, as only the regime F ^> 1 
was described for classical treatments of the system [7[ I15j. 

Figure 121 shows the time evolution of the moments (a + a^) and (n), obtained 
numerically, for a particular choice of parameters. For the relatively fast oscillator chosen 
(r = 0.6) the average resonator position shows a clear initial transient before decaying 
exponentially as if it were a damped harmonic oscillator. The average SET island charge, 
(n), relaxes on two time-scales: initially, the charge relaxes rapidly towards its steady 

§ We did check by numerical integration that for Gaussian initial states of the resonator, the master 
equation preserves positivity (even though it is not of the Lindblad form) and relaxes to a Gaussian 
steady state. 
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Figure 2. Evolution of the average resonator position, (a + aJ), and SET island 
charge, (n). Results from a numerical integration of the full master equation are 
compared with results calculated using the Caldeira-Leggett (CL) and Rotating Wave 
Approximation (RWA) master equations (described in sections 0] and respectively). 
Notice that the Caldeira-Leggett curves used 'slipped' initial conditions [21] chosen 
to fit the long-time behaviour of the numerical results. The parameters used are 
eV/hw = 6, g = 0.1/h (i.e. T = 0.6) and A = 0.2. 



state value over a time geV, it then undergoes weak oscillations (following those of the 
resonator) which are damped out slowly. 

Although the dynamics of the resonator is qualitatively the same for the whole range 
of parameters we have studied, there are important quantitative differences. In figures 
|3]andE]we plot the damping rate and frequency shift (extracted from fits to the average 
resonator position as a function of time) for different choices of eV, A and for a wide 
range of T values. Strikingly, we find that the data all fall very close to universal curves. 
When time is measured in units of the resonator period, the form of the universal curves 
for damping and frequency shift are determined by T and k = 2\ 2 huJo/eV = muj^x^/ eV 
which gives a measure of the SET-resonator coupling strength compared to to the bias 
voltage. As we shall see later on, the universal curves emerge naturally from a reduced 
master equation for the resonator. 

From figure EH we can see that the damping rate reaches a maximum when the 
resonator period matches the tunnelling rate, and drops to zero when the tunnelling rate 
is either much faster or much slower than the oscillator. We also note that the damping 
rate is proportional to V for T 1 and inversely proportional for V ^> 1, with the latter 
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Figure 3. Plot of the resonator damping rate as a function of the average tunnel 
rate, T. For simplicity, the junction resistances are assumed to be equal, gu = g$ = g, 
and we choose n' g — 1/2 (the degeneracy point). The voltages given for the various 
data sets are measured in units of tkuQ. Also shown are analytic curves which match 
expressions derived in sections 0] and [SJ 



matching the results obtained in previous classical treatments of the system E3 
while the former is consistent with results found using a quantum treatment in 
The frequency shift follows a different pattern, increasing monotonically with T until 
saturating at a constant value for r 3> 1. 

Physically, the frequency shift and damping of the resonator arise because the 
charge on the SET is a function of both the resonator position and (indirectly) velocity 
and the resonator in turn experiences a force proportional to the charge on the resonator. 
Thus the resonator-SET coupling results in forces on the resonator proportional to its 
position (frequency shift) and to its velocity (damping) [H]. The velocity dependent 
term is due to the fact that the SET takes a finite amount of time to respond to 
the resonator. If the SET responds instantly (which is equivalent to taking the limit 
r — > oo,) there is no damping. Conversely, if the SET takes an essentially infinite 
amount of time to respond to the resonator position (i.e. in the limit r — > 0), then the 
damping must also disappear. 

The second moments of the resonator can also be calculated from equation El and 
solved to find the steady state variances of the resonator. We find, 

m [ Su2) ST = mw 2 {5x 2 ) ST = eV(n) ST [l - (n) ST ] (21) 
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Figure 4. Plot of the relative frequency shift (squared) as a function of T. The 
junctions are chosen to have equal resistance go = 9s = g, and we choose n' g = 1/2 
(the degeneracy point). The voltages given for the various data sets are measured in 
units of hujQ. The analytic curve comes from equation (|24|l discussed in section 



where 

eV - 2E C (1 - 2n'\ 
<»>«* = 2eV(l > 5 ' < 22 » 



is the steady state value of the island occupation (given by equations ITSll2Tjj) . In the weak- 
coupling limit, n — > 0, equipartition is recovered and if we also take the high voltage 
limit, eV ^> Hujq, then|| we can identify an effective temperature for the resonator which 
matches previous classical results [Zj: ksTsET = eV{n)sT\\ — ( n )sr]- Note that these 
expressions for the steady state variances are not valid for arbitrarily small voltages 
(and hence do not violate the uncertainty relations) as they were derived assuming a 
bias voltage large enough to ensure that no back-tunnelling occurs. Specifically, equation 
(|2Tj) is valid for eV > 2(E ch + 2Huj ) and eV > 2(—E ch + 2ftw ), and for weak electro- 
mechanical coupling. 



4. Caldeira-Leggett Equation for the Resonator 

In section El we traced over the microscopic electronic degrees of freedom to obtain a 
master equation for a system consisting of the resonator and the macroscopic charge of 
the SET island. A natural further step is to trace over the remaining electronic degree 

|| Away from the high voltage limit, the value of the SET effective temperature will depend on the 
resonator frequency, as pointed out in 
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of freedom and hence obtain a reduced master equation for the resonator alone |15j . 
The advantage of this reduced master equation is that it gives a compact description 
of the resonator dynamics and makes clear the analogy between the effect of the SET 
electrons and a thermal bath. 

In order to trace over the macroscopic charge variable, we treat the SET as the bath 
and make many of the same kinds of assumptions as were used in deriving equation (JTfij) 
(details of the calculation are given in |Appendix B ) . In particular, we assume that 
the tunnelling rate, geV, is fast compared to the rate of change of the resonator's 
reduced density matrix in the interaction picture (Markov approximation). Whilst we 
cannot be completely sure when such an approximation is valid before carrying out 
the calculation, we can judge its validity by comparing the predictions of the resulting 
master equation with those of equation (fT6|) described in the previous section. It is 
also necessary to assume that the coupling is weak compared to the bias voltage (Born 
approximation), k -C 1, but this is not restrictive given the weak-coupling assumption 
used in the derivation of equation (flfijl. We also assume that the SET and resonator 
density matrices factorise. 

The reduced master equation for the resonator we arrive at by means of the 
Born-Markov approximations takes a very similar form to that of the Caldeira-Leggett 
equation which describes quantum Brownian motion [23 EH! I2Z] > 

p{t) = -i[H rjPr ( t )} - {p,Pr(t)}) ~ D[X, [x, Pr (t))) - l[x, \p, P r(t))) (23) 

where H r is the Hamiltonian of an harmonic oscillator with the renormalised frequency, 

_r_ 2 
TTr 2 



u R = u ( 1 - n^-^ ) (24) 



and the damping constant is^f 

kT 

Icl = uJo x - p • (25) 
The normal and anomalous diffusion constants fZG\ ITS] are 

D= m ^ eV (n) ST [l-(n) ST ] (26) 

and 

h f = -TTtW^M 1 - (»>ot], (27) 

(i + r 2 ) 

respectively. 

The steady state resonator variances are readily obtained from the reduced master 
equation and hence we obtain 

mu^Sx 2 ) ST = m(Su 2 ) ST — hf (28) 

m(5u 2 ) ST = eV(n) ST [l - (n) ST ] (29) 

mu 2 (5x 2 ) ST = eV(n) ST [l - (n) ST ] ! ^ • (30) 

1 i+r 2 

For comparison, the clamping and renormalized frequency found in previous classical calculations [3 
El El were 7 = kloq/T and lor = lo q {\ — n) 1 / 2 , respectively. 
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These values are somewhat different to those obtained before the macroscopic charge 
variable was traced out. In particular, the presence of the anomalous diffusion term 
means that the variances in the position and velocity are no longer direcly proportional 
to one another" 1- . Such differences must be interpreted as artifacts of our Born-Markov 
approximations. However, in the limit k — > these values converge to the same limits 
as before. 

Having compared the steady state properties of the resonator given by the reduced 
master equation with the 'true' values given by the master equation for the resonator 
and SET island charge, we now compare the dynamics. In figures |3] and 0] we plot 
the analytic Caldeira-Leggett values of frequency shift and damping [equations (J23|) 
and (|25p] together with values extracted from numerical integrations of the full master 
equation. The analytic results match up very well with the numerics in the weak electro- 
mechanical coupling limit. Furthermore, this good agreement extends over the whole 
range of T showing that our Markov approximation remains valid even for r< 1. 

Although the reduced master equation successfully describes the long-time 
dynamics of the resonator, it fails to capture certain transient features which are present 
in the numerics over time-scales of order the electron tunnelling time, geV. In figure El 
we compare the time evolution of the moments calculated from equation to those 
obtained from the full master equation. In order to match the long time behaviour of the 
numerics we have had to 'slip' the initial conditions for the reduced master equation [24J. 
The failure of a reduced master equation to capture transient features in the motion of 
the system is often a consequence of Markov approximations used in their derivation [24J, 
and it seems very likely that the same has happened here. 

An additional factor that may introduce further complications is the presence of 
an external bath, which has been neglected in this treatment. While naively we might 
expect that a thermal environment could simply be incorporated into the reduced master 
equation of the resonator (equation EH} by the addition of further Caldeira-Leggett type 
terms the correlation time of the environment provides an additional timescale 
which might affect the validity of the assumptions used in deriving the reduced master 
equation in the first place [2E1- A more detailed treatment would therefore required, for 
instance including the environment at a more microscopic level as a bath of harmonic 
oscillators in the original Hamiltonian of the system. 

5. Rotating Wave Approximation 

An alternative way of deriving a reduced master equation for the resonator alone is 
to make use of a rotating wave approximation (2TJ, |22j (RWA). In the RWA, terms 
with an e luJot dependence are assumed to oscillate very rapidly compared to the rate of 
change of pi and hence effectively average to zero. As with the validity of the Markov 

+ The presence of an anomalous diffusion term in the reduced master equation of a resonator coupled 
to a point contact was reported in [T^]. However, in that case there is no analogue of the resonator-SET 
charge master equations which can be used to determine the 'true' dynamics of the resonator. 
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approximation used in deriving the Caldeira-Leggett type master equation, we cannot 
be sure in advance whether this condition will be satisfied, but we can verify that it is 
a posteori from the rate of change of pi in the resulting master equation. 

Returning to equation (j!5|) . implementing the RWA means dropping any term where 
(ui m — io n ) 7^ 0. Expanding up to order A 2 , the resulting master equation is 

Pj = gsEsT [cr~, cr + {l — \ 2 (a)a + aa^)}] pi 

+ g s {(E s - hu^X 2 ? [aV - , aa + ] + (E s + frw^X 2 ? [aa~, aV + ] } p, 
+ gr>Ei£>T [o~ + , cr~ {l — \ 2 (o)a + aa^)}] pi 

+ g D { (E D - ftwa)\ 2 F [a ] a + , aa~] + (E D + Huj )X 2 J r [aa + , aV~] } pi. (31) 

This master equation is closely related to that found in [Oj where a master equation for 
a closely related system of a quantum dot gated by a nanomechanical resonator was 
derived. 

When we make the simplifying assumption g$ = go — g, the dynamics of the 
resonator and SET island charge decouple. The equation of motion for the average 
charge then takes the simple form, 

j t (n) = -g [E s - hio^X 2 ] (n) +g[E D - fru \ 2 ] (1 - (n)). (32) 

The behaviour of (n) according to this equation is compared with the 'true' behaviour 
given by equations (fTHj) -(|2T? j) in Figure El It is clear that although the fixed point value 
of (n) remains the same, the dynamics is different as it no longer displays oscillations 
correlated to the resonator motion. 

Assuming that the density matrix factorises into charge and resonator parts, the 
evolution of the resonator alone is given by a reduced master equation which has the 
Lindblad form 

P r ,i = 9^ 2 {(n)(E s - E D ) + (E D - hwo)} F[a\ a]p rJ 

+ g\ 2 {(n)(E s - E D ) + (E D + hu )} F[a, <f]p r j, (33) 

where 

p r ,i = (N \pi\N ) + (JVi|joj|JVi>. (34) 

The average charge appears in equation ()33|) as a parameter though it is not necessary 
to assume that (n) has reached its fixed point value. Instead, we can use the solution 
of equation (J3~2"|) to obtain a reduced master equation for the resonator capable of 
capturing transient behaviour arising from different initial conditions, albeit one with 
time dependent coefficients. 

The decoupling between the occupancy of the SET island and the resonator motion 
occurs because we are working in the limit where the resonator oscillations are much 
faster than the time-scale for electron tunnelling events and hence the SET only 
experiences the average of the resonator position. It should be noted, however, that 
this decoupling only occurs in the strict RWA limit, i.e. r — >• 0: for any non-zero values 
of T there is a residual correlation between the resonator motion and the average charge. 
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The oscillations in (n) [obtained from equations (|18| ) - (j2(Jj) ] shown in Figure 121 for which 
we set T = 0.6, clearly follow those of the average resonator position, but this behaviour 
disappears completely in the RWA limit. 

Equations of motion for the moments of the resonator are readily obtained from 
equation or from equation (JHH), 

^a) + a) = -iu (a-a ] ) - g\ 2 huj (a ] + a) +2g\ 3 hw (n) + 2\^-(n) (35) 

(JjL (JjL 

i— (at — a) = — ujq(o) + a) + 2Aco>o(n) — ig\ 2 frwo(a^ — a). (36) 

In order to capture the long-time dynamics of the resonator as simply as possible, it is 
sufficient to assume that the average island charge has reached its steady state values. 
We can then obtain the following equation of motion for the resonator position about 
its fixed point value, 

(x) = —u 2 (x) — 2g\ 2 fkjQ(x) , (37) 

where we have discarded terms 0(X 3 ) and higher. From this equations it is clear that 
the resonator is damped at a uniform rate ^rw = 2g\ 2 huJo = kTuq, but there is no shift 
in its frequency. The steady state values of the resonator variances take the form, 

mu^Sx 2 ) ST = eV(n) ST [l — (n) ST \ + \ 2 m<J^x 2 q (38) 
m(8u 2 ) ST = (eV — A\ 2 muj 2 x 2 q ) (n) ST [1 — (n) ST \ + X 2 muj1x 2 q . (39) 

As with the variances obtained from the Caldeira-Leggett master equation, these differ 
slightly from the 'true' values [equation but they do approach the same weak- 

coupling limit (A — > 0). 

In figure 121 we compare the RWA evolution of the system moments with the full 
case. We see that the agreement is very close, including the short-time transient, for 
the value T = 0.6 used and in fact it becomes exact in the limit r — > 0. The RWA 
values for the damping and temperature, (and the absence of a frequency shift) are 
consistent with those calculated from equations (fTKJ) and (fTT)J) in the limit r < 1, as 
can be clearly seen from figures HU and EJ This implies that the correct criteria for the 
validity of the rotating wave approximation in this system is T 1. This is in contrast 
to the case of just a harmonic oscillator coupled to a bath j^I], where the criteria for 
the validity of the RWA is simply that the resonator frequency should be much larger 
than the damping rate due to the bath. It is clear that the presence of the SET in our 
system introduces an additional time-scale over which pi evolves and hence the period 
of the resonator oscillations must be much shorter than both the damping time and the 
electron tunnelling time for the RWA approximation to be valid. 

It is worth noting that equation EH] seems to suggest that the damping is caused by 
terms like T acr~] pi which add or remove an energy quantum from the resonator 

during a tunnelling event (for discussions of energy exchange in similar master equations 
see jSlEHl)- Although this mechanism seems intrinsically quantum mechanical, involving 
the exchange of energy in units of fujQ, the RWA master equation was derived from an 
equation describing essentially classical behaviour. However, we must remember that the 
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equation we have derived is an unconditional master equation. Even though the terms 
in equation EU describe the transfer of individual quanta, the fact we are considering 
an ensemble of systems means that the density matrix changes continuously. A full 
description of damping that proceeds by transfer of discrete units of huj Q (as equation EH1 
seems to imply) would require quantum trajectory methods |22j- 

6. Conclusions and Discussion 

Starting from a microscopic formulation, we have derived a quantum master equation 
for the SET-resonator system which describes the coupled dynamics of the SET 
island charge and the resonator in the limit of weak electro-mechanical coupling and 
intermediate voltages (i.e. eV > ftw ). By investigating the dynamics of the charge and 
resonator moments we find that the resonator is damped, undergoes a shift in frequency 
and eventually reaches a thermal-like steady state. These results agree with those found 
in a classical description of the same problem [TJEi], are more general. In particular, 
we find that the resonator is damped by the SET even when the electron tunnelling 
time is of order the resonator period or longer. Furthermore, maximum damping of the 
resonator by the SET electrons occurs when the electron tunnelling time matches the 
resonator period. 

Starting with the resonator-island charge master equation, we have derived reduced 
master equations for the resonator alone in two different ways. The first reduced master 
equation is very similar to the Caldeira-Leggett master equation and is obtained by 
tracing over the island charge, treating it like an external bath (i.e. making Born- 
Markov approximations). The second type of reduced master equation is found by 
making the rotating wave approximation so that the charge dynamics decouples from 
the resonator motion. Apart from an initial transient, the Caldeira-Leggett type master 
equation describes the long-time dynamics of the resonator faithfully over the whole 
range of relative electrical and mechanical time scales. In contrast, the master equation 
obtained via the rotating wave approximation only captures the resonator dynamics 
when its period is much shorter than the electron tunnelling time, though unlike the 
Caldeira-Leggett master equation it captures the transient motion and is of the Lindblad 
form so is guaranteed to preserve the positivity of the density matrix. 

The derivation of a reduced master equation for the resonator in the SET-resonator 
system makes an interesting comparison with better known dissipative systems as it 
proceeds via a two stage process: tracing over the microscopic electronic levels gives 
rise to a master equation which still contains the SET island charge, and further 
approximations are needed to obtain a description of the resonator alone. This is 
in contrast to other NEMS such as the resonator-point contact system [HI EJ ITH] 
where tracing over the microscopic electronic states leads directly to a Caldeira-Leggett 
type master equation for the resonator. However, the intermediate stage for the 
SET-resonator system (consisting of a master equation for the resonator and SET 
island charge) leads rather fortuitously to closed sets of equations of motion for the 



Quantum master equation descriptions of a nanomechanical resonator 



18 



resonator and charge moments. Thus it is possible to characterize the dynamics of the 
resonator using simple numerical and analytical techniques before making the additional 
approximations required to obtain the reduced master equations of the resonator. Thus 
we are able to show that the anomalous diffusion term which arises in the Caldeira- 
Leggett type reduced master equation is not a feature of the 'true' resonator dynamics, 
but rather an artifact of the approximations we make. 

The quantum master equations we derive here for the SET-resonator system will 
provide useful tools for further investigations into the quantum dynamics of this system. 
In particular, we can use them as a starting point for investigations of how an individual 
SET and resonator (rather than an ensemble of such systems) evolves. However, it will 
also be interesting to try to extend the ensemble-averaged master equations derived 
here to obtain descriptions of the low-bias and strong-coupling regimes where the close 
connection between quantum and classical treatments found here is unlikely to persist. 
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Appendix A. Alternative derivation of the resonator and island charge 
master equation 

In the body of the paper we rewrote equation to get a master equation in a simple 
form (equation [THJ) • Although the derivation of equation [T3] required an expansion in A, 
this simple master equation is in fact more general. Here we carry out a slightly different 
derivation that does not require A < 1, and is therefore valid in the fully classical limit 
where both the coupling \ an d the bias voltage V are much larger than ftw . Instead, 
we make a different weak coupling approximation and assume that the coupling is small 
compared to the bias voltage, x <S eV. In addition, we assume that the resonator moves 
slowly on the relaxation time of the electrons in the leads. 
We start from the Born-Markov master equation, 



PAt) =-^J dt '^ [Bi®, [#/(0,p/(*)s]] 





where the interaction picture Hamiltonian (equation EJ), is 

Hj = J^Tu [4(t)c t (t)a-(t)e^ at ^ ot - ae ^ + c^cl^^e-^^^'^] . (A.l) 
i 

Using the Baker- Campbell-Hausdorff theorem [21] to rearrange the terms under the 
integral so that the time dependence in the exponentials is expressed in terms of r = t—t', 
we then assume that the resonator evolves slowly on the time-scale of the bath correlation 
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functions (i.e the relaxation time of the electrons in the leads) and hence make the 
approximation e luJ ° T ~ 1 + iuj t. Then we perform the integration over r to obtain 



Pi® = ^El T ^4^ \s'Ht)e' X ^- a \S'(t) 



Li 



x5 (a - e, + E ch - A^ (a f + a) - \ 2 hu ) e x ^~ a ^ pi(t)E 
+ ^Tr^|^| 2 4s a jr [5'(t)e A ( at - a ),5't(t)e- A ( at - a ) 

x5 (ej - ei + E ch - Xtko Q (a) + a) - X 2 fku )] pj(t)E. (A.2) 

where S'(t) = a~ e~ l ( Ech / ti ~ >J * > ( a +a)-A 2 w )* _ Although this appears unwieldy, converting 
back from the interaction and canonically transformed picture using the definition, 
p(t) = e A ( at ~ a )e~* Hs *p/(t)e +iiTs *e~ A ( at ~ a ), the above equation reduces to a simple form. 
After tracing over the bath, performing the integrals over the Fermi functions (following 
the same approach as in section [2} , and assuming that the electro-mechanical coupling 
is sufficiently weak to ensure that no back-tunnelling occurs, we regain equation ED 

pit) = - l - [H s , pit)} + g s T [a~(E s - \hw Q {a) + a) + \ 2 hw ), a + ] p(t) 

+ g D T [a + (E D + \huJo(a^ + a) - X 2 hu ), a~] p(t). 

Note that in this case the condition on the coupling required to ensure that no back- 
tunnelling occurs is A <§; eV/hcj (i.e \ ^ eV)- 

Writing the charge- diagonal matrix elements (N \p\No) and (Ni\p\Ni) of equation 
(HHJ) in Wigner function form [2JJ results in a pair of coupled equations equivalent to 
the classical master equations proposed in ref. [7j 
/ q q\ 

Woa(x,u) = \ ulx— -u-q^J W 00 (x,u) + g s [E s - mu 2 x x + X 2 huj ] W n (x,u) 

— go [Ed + muj 2 XQX — X 2 hwo\ Wqq(x, u) (A. 3) 

/ d d\ 

W n (x,u) = {uJq(x - xq)— - u -q^J W u (x,u) - g s [E s - mu 2 x x + \ 2 fixJo] W u (x,u) 

+ 9d [Ed + mLU 2 x x - X 2 fku ] W 00 (x, u). (A.4) 
Appendix B. Caldeira-Leggett equation for Resonator 

It has been shown 15] that when the resonator coupling is small compared to the gate 
voltage (x <C eV), and the separation of timescales is large (r ^> 1), a further Born- 
Markov approximation can be performed on the classical master equation jjj, to get a 
Fokker-Planck equation for the resonator alone. We can perform a closely analogous 
procedure, starting with equation [TBI to get a master equation for the resonator alone, 
but we do not make the assumption r ^> 1. 

First, we note that the diagonal elements of the SET p 00 = (N Q \p\N ) and 
pn = (Ni\p\Ni) decouple from the off-diagonal terms. Writing the density matrix 
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as a vector p = (poo>Pn), we can re- write equation ITol as. 

i r . I (a f + a) 




+ »("I D -l-J'+^f" 1 ) {(»'+«). rf. (B-D 

where if r = ftco>o a ^, the curly braces, {•,•}, indicate an anti-commutator and we have 
shifted the origin of the oscillator position by the steady state value (a) + a) ST so that 
the new fixed point will be zero. The electrostatic energy differences associated with 
tunnelling are given by 

eV 

E' D = E D - {\ 2 hujQ + \hwo{a) + a) ST ) = — E c (l - n' g ) + Xhw (a^ + a) ST (B.2) 

eV 

E' s = E s + (\ 2 hw + \huj {a ] + a) ST ) = — + E c (l - n' g ) - \huj (a^ + a) ST . (B.3) 

To obtain equation [THl we traced over the electron levels to get an effective 
master equation for the resonator and the charge on the SET island. Now, we wish 
to trace over the remaining charge degree of freedom to leave a master equation for 
just the resonator. We switch to an interaction picture representation defined by 

p(f) = e 9H S ETt e -iH r t/h p ^ e iH r t/h^ where 



H 



SET 



~ e 'd e 's 
e 'd -E's 



is the non-Hermitian 'Hamiltonian' that describes the non-unitary evolution of the SET 
in the absence of interaction with the resonator. We define the SET-only operators as 

Bl = 9 -^[ * I (B.4) 




b-2 /^,( x 2 /M - a( ) ). <b^ 

In this interaction picture, the equation of motion for the density matrix is 

pi{t) = - l -B 2 (t)[{a) + a)(t), Pl (t)} + Si(t){(at + a)(t), Pl (t)}. (B.6) 

Following a very similar procedure to that used to obtain equation (jA.l|) . we perform 
a Born-Markov approximation ^H] and hence obtain an equation of motion for the 
resonator alone, 

t 

p'n(t)= J TrB 2 (t)B 2 (t')PsET(0)[(ai + a)(t),[(ai + a)(t'),pi(t)]]dt' 

- % - I M 2 (t)S 1 (OPs Br (0)[(ot + o)(*),{(aU «)(/'). ,>,(t)}} <lt' 
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TrB 1 (t)B 2 (t')P SET (0){(a^ + a)(t),[(a t + a)(t'), Pl (t)}} dt' 

o 

t 

+ j TrB 1 (t)B 1 (t')P SET (0){(a^ + a)(t),{(a^ + a)(t'),p I (t)}}df (B.7) 
o 

where Pset is a vector that describes the occupation of the SET. The Born 
approximation is based on the assumption of weak electro-mechanical coupling which 
from equation (jB.lj) implies huioX 2 <C eV (i.e. k < 1). The Markov approximation 
assumes that the SET correlation functions decay very rapidly on the time-scale of the 
rate of change of the interaction picture density matrix [2*2*] . 

We now need to evaluate the SET correlation functions. The only non-zero terms 

are, 

Tr[B 2 (t) Bl (t')P SET (0)] = ~ ^f^ e-^vit-t') (Rg) 

Tr[B 2 (t)B 2 (t')P SET (0)} = (hLOofX'^^e-^-^, (B.9) 

and, in line with our Markov approximation, we discard terms that decay like e~ geVt * . 
Noting that (a* + a)(t) = (aJ + a) cosfwot) + i(aJ — a) sin(uoi) /ujq, we can insert the 
correlation functions into equation IB .71 and do the integrals over t' exactly. Converting 
back to the Schrodinger picture, 

p r (t)= - l -[H r ,p r {t)) 

1 f X 2tf Uo E D E 's\ [(^ + a),[T(a^ + a)-t(a^-a),p r (t)]] 



h 2 V (eV) 2 J T 2 + 1 

i g\ 2 h 2 uj [(a f + a), {r(a t + a) - i(at - a),p r (t)}] 

+ h—2 pTI ' (R10) 

we find a master equation that can be easily rearranged into a more familiar form, which 
is very close to the Caldeira-Leggett equation, 

P® = - \[Hr,Pr{t)] ~ \ ( ^° i+ r r2 ) [X 2 ,pr(t)] ~ % ^f[x, { P , Pr (t)}} 

-D[x,[x,p r (t)}]-£[x,\p,p r (t)}}. (B.ll) 

The second and first terms can be combined, leading to a renormalization of the 
resonator frequency, 



U R = U Q [1- K 



r 



2 



T + r 2 

and the damping, diffusion and anomalous diffusion constant are given by 
2\ 2 gftw kT 



* It is therefore no surprise that the equation we eventually obtain shows transient deviations from the 
'true' resonator behaviour on times ~ geV. 
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D 



m-fcLE' D E' e 



s 



eVH 2 

hf = 



E' d E' s k 



eV(l + Y 2 y 

respectively. These values are consistent with previous (classical) calculations in the 
limit r > 1 [2EEH]. 
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